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Abstract. Exotic dark matter together with the vacuum energy (cosmological con- 
stant) seem to dominate in the flat Universe. Thus direct dark matter detection is 
central to particle physics and cosmology. Supersymmetry provides a natural dark 
matter candidate, the lightest supersymmetric particle (LSP). Furthermore from the 
knowledge of the density and velocity distribution of the LSP, the quark substructure of 
the nucleon and the nuclear structure (form factor and/or spin response function), one 
is able to evaluate the event rate for LSP-nucleus elastic scattering. The thus obtained 
event rates are, however, very low. So it is imperative to exploit the two signatures of 
the reaction, namely the modulation effect, i.e. the dependence of the event rate on the 
Earth's motion, and the directional asymmetry, i.e. the dependence of the rate on the 
the relative angle between the direction of the recoiling nucleus and the sun's velocity. 
These two signatures are studied in this paper employing various velocity distributions 
and a supersymmetric model with universal boundary conditions at large tanj3. 

1 Introduction 

The combined MAXIMA-1 BOOMERANG and COBE/DMR Cosmic Microwave 
Background (CMB) observations imply that the Universe is flat pi, Q — 1.11 ± 
0.07, while the baryonic component is very small H^h 2 — 0.032 Q 0u g 9 . Further- 
more exotic dark matter has become necessary in order to close the Universe. In 
fact about a decade ago the COBE data |2| suggested that CDM (Cold Dark 
Matter) dominates the Universe, Qcdm being at least 60% Subsequent ev- 
idence from two different teams, the High-z Supernova Search Team [|| and the 
Supernova Cosmology Project || : Q changed this view suggesting that the 
Universe may be dominated by the cosmological constant A or dark energy. In 
other words one roughly finds a baryonic component fls = 0.1 along with the 
exotic components ficDM — 0.3 and Qa = 0.6. In a more detailed analysis by 
Einasto § one finds h = 0.65 ± 0.07, fl h = 0.05 ± 0.02, Q C dm = 0.36 ± 0.05, 
^hdm < 0.05 and Qa — 0.65 ± 0.05. Since the non exotic component cannot 
exceed 40% of the CDM [§, there is room for tnc exotic WIMP's (Weakly 
Interacting Massive Particles). In fact the DAMA experiment |9j has claimed 
the observation of one signal in direct detection of a WIMP, which with better 
statistics has subsequently been interpreted as a modulation signal []Tof . 

In the most favored scenario of supersymmetry the LSP can be simply de- 
scribed as a Majorana fermion, a linear combination of the neutral components 
of the gauginos and Higgsinos ]Tl| , p^| , p^[ . 
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2 An Overview of Direct Detection - The Allowed SUSY 
Parameter Space. 

Since this particle is expected to be very massive, m x > 30GeU, and extremely 
non relativistic with average kinetic energy T < lOOA'eU, it can be directly 
detected [|l6[0 mainly via the recoiling of a nucleus (A,Z) in elastic scattering. 
In order to compute the event rate one needs the following ingredients: 

1) An effective Lagrangian at the elementary particle (quark) level obtained 
in the framework of supersymmetry as described , e.g., in Refs. [114 f 



2) A procedure in going from the quark to the nucleon level, i.e. a quark model 
for the nucleon. The results depend crucially on the content of the nucleon in 
quarks other than u and d. This is particularly true for the scalar couplings as 
well as the isoscalar axial coupling fH^ — [El] . 

3) Compute the relevant nuclear matrix elements using as reliable as 
possible many body nuclear wave functions. The situation is a bit simpler in the 
case of the scalar coupling, in which case one only needs the nuclear form factor. 

Since the obtained rates are very low, one would like to be able to exploit 
the modulation of the event rates due to the earth's revolution around the sun 
[ p5|j2^| — [p8| . To this end one adopts a folding procedure assuming some distribu- 
tion p6pq of velocities for the LSP. One also would like to know the directional 
rates, by observing the nucleus in a certain direction, which correlate with the 
motion of the sun around the center of the galaxy and the motion of the Earth 

The calculation of this cross section has become pretty standard. One starts 
with representative input in the restricted SUSY parameter space as described 
in the literature Jl2|,[l4| (see also Arnowitt and Dutta which will appear in 
these proceedings) . We will adopt a phenomenological procedure taking universal 
soft SUSY breaking terms at Mqut, i-e., a common mass for all scalar fields too, 
a common gaugino mass M1/2 and a common trilinear scalar coupling Aq, which 
we put equal to zero (we will discuss later the influence of non-zero Aq's). Our 
effective theory below Mqut then depends on the parameters O]: 

mo, Mi/2, Mo, olg, Mqut, h t , , h, , h T , tan/3 , 

where etc — 9%/^ (dG being the GUT gauge coupling constant) and h t ,hb, h T 
are respectively the top, bottom and tau Yukawa coupling constants at Mqut- 
The values of «g and Mqut are obtained as described in Ref . Q . For a specified 
value of tan/3 at Ms, we determine h t at Mqut by fixing the top quark mass 
at the center of its experimental range, m t (m t ) — 166GeV. The value of h T 
at Mqut is fixed by using the running tau lepton mass at mz, m T (mz) = 
1.746GeV. The value of hb at Mqut used is such that: 

^MmzjfE = 2-90 ± 0.14 GeV. 



after including the SUSY threshold correction. The SUSY parameter space is 
subject to the following constraints: 
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1. ) The LSP relic abundance will satisfy the cosmological constrain: 

0.09 < f2 L sph 2 < 0.22 (1) 

2. ) The Higgs bound obtained from recent CDF |cj and LEP2 |l[ , i.e. m h > 113 

3. ) We will limit ourselves to LSP-nucleon cross sections for the scalar coupling, 
which gives detectable rates 

4 x 1CT 7 pb < o n s ^l e a ° n < 2 x 1CT 5 pb (2) 

We should remember that the event rate does not depend only on the nucleon 
cross section, but on other parameters also, mainly on the LSP mass and the 
nucleus used in target. The condition on the nucleon cross section imposes severe 
constraints on the acceptable parameter space. In particular in our model it 
restricts tanj3 to values tan/3 ~ 50. We will not elaborate further on this point, 



since related work has already appeared elsewhere 112 32 



3 Expressions for the Differential Cross Section . 



The effective Lagrangian describing the LSP-nucleus cross section can be cast 
in the form [M 



Gp 
V2 



{(Xi7 A 75Xi)^a + (xiXi)^} 



(3) 



where 



Jx=N lx {f v + f v T 3 + f Al5 + f 1 Al5 T 3 )N , J = N(f° + flr 3 )N (4) 

We have neglected the uninteresting pseudoscalar and tensor currents. Note 
that, due to the Majorana nature of the LSP, Xi7 A Xi — (identically). 

With the above ingredients the differential cross section can be cast in the 
form 0,||2G 



da(u, v 
£s = 00 ( 



+ S V ^) F\u) + £ spin F u (u)} 



X A 2Z 2 s . rr \->-! 
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s(u) = [(^n (o))^-pfi + 2^(0)^ (o)^ + ^i(o)) 2 - 



Fn(u) f\ 



Fu(u) 



(5) 
(0) 
(7) 

(8) 
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Zv = < Y c Cv (9) 



(fi r (A)/ /j, r (N)) 2 a2 n fiA-2Z 2 v 2 1 2r7+l ( 2u 

= n A U - 77f ^ J LC—J i 1 



(1 + 4)2 V ft A c (2 Mr 6)2(l + ? /) 2 (" 2 ) J 

J V 

(10) 

°"p x° = P r °ton cross-section, i = S, spin, V given by: 

a px° = CT ° (fs) 2 ) 2 (scalar) , (the isovector scalar is negligible, i.e. er^ = 

4) ' 

<£5 = 3 (/I + /i) 2 (^f) 2 (spin) , < x0 = a (ft+ftf (^f) 2 (vector) 
where is the nucleon mass, 77 = m x /m]^A, and jtx r (-A) is the LSP-nucleus 
reduced mass, fi r (N) is the LSP-nucleon reduced mass and 

co = ^(Gftoat) 2 ~ 0.77 x 10" 38 cm 2 (11) 

Q = Q u , Q = . 1 u9 = 4.1 x 10 4 A- 4 / 3 KeF (12) 

where Q is the energy transfer to the nucleus and F(u) is the nuclear form factor. 

In the present paper we will concentrate on the coherent mode. For a discus- 
sion of the spin contribution, expected to be important in the case of the light 
nuclei, has been reviewed elsewhere 24,32|. 



4 Expressions for the Rates. 

The non-directional event rate is given by: 

n r, dN p(0) m 

H = Knon-dir = —rr = < cr(u,V)\v\ (13) 

at m x Aitin 

Where p(Q) — O.SGeV/cm 3, is the LSP density in our vicinity and m is the 
detector mass The differential non-directional rate can be written as 

dR = dR non ^ dlr = m da(u,v)\v\ (14) 

m x Amff 

where da(u,v) was given above. 

The directional differential rate [|ll] , [^8| in the direction e is given by : 

p(0\ 7TL 1 

dRdir = : v.eHtv.e) — da(u,v) (15) 

m x Aitin 2tt 

where H the Heaviside step function. The factor of 1/2-7T is introduced, since the 
differential cross section of the last equation is the same with that entering the 
non-directional rate, i.e. after an integration over the azimuthal angle around 
the nuclear momentum has been performed. In other words, crudely speaking, 
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l/(27r) is the suppression factor we expect in the directional rate compared to 
the usual one. The precise suppression factor depends, of course, on the direction 
of observation. The mean value of the non-directional event rate of Eq. ([l4j), is 
obtained by convoluting the above expressions with the LSP velocity distribution 
J(v,ve) with respect to the Earth, i.e. is given by: 

/ dR\ p(0) m f „, . . ,da(u,v) . „, 

<^>=^i^y^^H^^ (i6) 

The above expression can be more conveniently written as 

/dRv Ko) m ^- ds ds r m /( ViVE) ^l d s v (17) 

\du/ m x Am N v v n du' v du' J yftv 1 ) du y J 

After performing the needed integrations over the velocity distribution, to 
first order in the Earth's velocity, and over the energy transfer u the last expres- 
sion takes the form 

R = Rt[l + h(a, Q mm )cosa)] (18) 

where a is the phase of the Earth (a = around June 2nd) and Qmin is the 
energy transfer cutoff imposed by the detector. In the above expressions R is the 



rate obtained in the conventional approach |16| by neglecting the folding with 
the LSP velocity and the momentum transfer dependence of the differential cross 
section, i.e. by 

m x ArriN c z 

where — S,V,spin contain all the parameters of the SUSY models. The 
modulation is described by the parameter h . 

The total directional event rates can be obtained in a similar fashion by 



suitably modifying Eq. (16) 



f f{V:VE) !LH^^l d ^ (20) 

\ du I dir m x AmN J 27T du 

The integration of the above equation is difficult. So we find it convenient to go 
to a coordinate system in which the polar axis is in the direction of of observation 
e, which in the above coordinate system is specified by the polar angle and 
the azymouthal angle <P. In this new coordinate system polar angle specifying 
the velocity vector is simply restricted to be < 6 < tt, while the azymouthal 
angle <p is unrestricted. Thus the unit vectors along the new coordinate axes, 
X, Y, Z, are expressed in terms of the old ones as follows: 

Z = sin cos &x + sin sin <I>y + cos 0z. (21) 

X = cos cos <Px + cos sin <Py — sin 0x. (22) 
Y = - sin<Px + cos<Py. (23) 
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Thus the LSP velocity is expressed in the new coordinate system as: 

v x = sin cos <Pvx + sin sin <Pvy + cos 0vz, (24) 



•v 



= cos cos <Pvx + cos sin <Pvy — sin 0vz , (25) 
v z = — sin#fx + cos^uy, (26) 



with ux = v sin cos <f>, = t> sin 8 sin (f>, vz — v cos 8. It is thus straightforward 
to go to polar coordinates in velocity space and get: 

/dR\ = p(0)^_ f Vm v s dv f 1 ^ f 27r d J(0 1 ^v^,^v E )da(u,v) 



\ du I dir m x Anix J J J 2ir du 

(27) 

Now the orientation parameters and <P appear explicitly and not implicitly 
via the limits of integration. The function / can be obtained from the velocity 
distribution, but it will not be explicitly shown here. Thus we obtain: 

Rdir = R[(tdir/2n) [1 + (hi - h 2 )cosa) + h 3 sina] (28) 

where the quantity tdir provides the non modulated amplitude, while /ii,/i2 
and /13 describe the modulation. They are functions of and 4> as well as the 
parameters a and Q m in- The effect of folding with LSP velocity on the total rate 
is taken into account via the quantity tdir-, which depends on the LSP mass. All 
other SUSY parameters have been absorbed in R. In the special case previously 
studied, i.e along the coordinate axes, we find that: a) in the direction of the 
sun's motion hi = h 3 = 0, b) along the radial direction (y axes) /13 = and c) 
in the vertical to the galaxy h 2 = 0. Instead of tdir itself it is more convenient 
to present the reduction factor of the non modulated directional rate compared 
to the usual non-directional one, i.e. 

fred = 5^ = W(2?r t) = K/(27r) (29) 

It turns out that the parameter k, being the ratio of two rates, is less dependent 
on these parameters. Another quantity, which may be of experimental interest 
is the asymmetry \R d i r (-) - Rdir (+)|/ {Rdir (-) + Rdir{+))- This is completely 
independent of all other parameters except the LSP mass and the velocity dis- 
tribution. The directional rates exhibit interesting pattern of modulation. Given 
the functions hi(a, Q m in), I — 1,2,3, one can plot the the expression in Eqs ( p"8| ) 
and as a function of the phase of the earth a. 



5 Results and Discussion 



The three basic ingredients of our calculation were the input SUSY parameters, 
a quark model for the nucleon and the velocity distribution combined with the 
structure of the nuclei involved. We will focus our our discussion on the coherent 
scattering and present results for the popular target 127 1. We have utilized two 
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nuclcon models indicated by B and C, for their description see our previous 
work p2[ , which take into account the presence of heavy quarks in the nucleon. 
We also considered the effects on the rates of the energy cut off imposed by the 
detector, by considering two typical cases Q m in = 0, 10 KeV. The thus obtained 
results for the non modulated non directional event rates, Rt, in the case of the 
symmetric isothermal model for a typical SUSY parameter choice |l2| ] are shown 
in Fig. El The two relative parameters, i.e. the quantities t and h, are shown in 



140 




m LSP (GeV) 



Fig. 1. The Total detection rate per (kg — target)yr vs the LSP mass in GeV for a 
typical solution in our parameter space in the case of 127 1 corresponding to model B 
(thick line) and Model C (fine line). For the definitions see text. 



Figs H and H respectively in the case of isothermal models. The case of non 




jSP mass— > 



Fig. 2. The dependence of the quantity t on the LSP mass for the symmetric case 
(A = 0) on the left as well as for the maximum axial asymmetry (A = 1) on the right 
in the case of the target 127 1. For orientation purposes two detection cutoff energies 
are exhibited, Qmin = (thick solid line) and Qmin = 10 keV (thin solid line). As 
expected t decreases as the cutoff energy and/or the LSP mass increase. We see that 
the asymmetry parameter A has little effect on the non modulated rate. 
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Fig. 3. The same as in Fig. g| for the modulated amplitude. We see that the modulation 
is small and decreases with the LSP mass. It even changes sign for large LSP mass. 
The introduction of a cutoff Qmin increases the modulation (at the expense of the total 
number of counts). It also increases with the asymmetry parameter A. 



isothermal models, e.g. caustic rings, is more complicated [M and it will only 
briefly be discussed here. 

In the case of the directional rates we calculated the reduction factors and 
the asymmetry parameters as well as the modulation amplitudes as functions 
of the direction of observation, focusing our attention along the three axes, i.e 
along +z, —z, +y, —y, +x and — x p7j . Since f rec i is the ratio of two parameters, 
its dependence on Qmin and the LSP mass is mild. So we present results for 
Qmin = and give an average as a function of the LSP mass (see Table 

As we have seen the modulation can be described in terms of the parameters 
hi, i = 1, 2, 3 (see Eq. (p8|)). If the observation is done in the direction opposite 
to the sun's direction of motion, the modulation amplitude hi behaves in the 
same way as the non directional one, namely h. It is more instructive to consider 
directions of observation in the plane perpendicular to the sun's direction of 
motion (0 = 7r/2) even though the non modulated rate is quite reduced in 
this direction. Along the — y direction (<£ = (3/2)7r) the modulation amplitude 
hi — hi is constant, —0.20 and —0.30 for A = 0, 1 respectively. In other words it 
is large and leads to a maximum rate in December. Along the +y direction the 
modulation is exhibited in Fig. 0. 



6 Conclusions 

In the present paper we have discussed the parameters, which describe the event 
rates for direct detection of SUSY dark matter. Only in a small segment of the 
allowed parameter space the rates are above the present experimental goals. We 
thus looked for characteristic experimental signatures for background reduction, 
i.e. a) Correlation of the event rates with the motion of the Earth (modulation 
effect) and b) the directional rates (their correlation both with the velocity of 
the sun and that of the Earth.) 

A typical graph for the total non modulated rate is shown Fig. [l]. The relative 
parameters t and h in the case of non directional experiments are exhibited in 
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Fig. 4. On the left figure one sees the quantity hi — hi in the direction +y for A = 
(thick line) and A = 1 thin line. In the —y direction this quantity is constant and 
negative, —0.20 and —0.30 for A = and 1 respectively. As a result the modulation effect 
is opposite (minimum in June the 3nd). On the right figure we show the modulation 
amplitude in the direction +x, which is essentially /13, since |/ii| << | /i3 1 - Thus in this 
case the maximum occurs around September the 3nd and the minimum 6 months later. 
The opposite is true in the —x direction. 



Table 1. The ratio k of the non modulated directional rate along the three directions 
to the non-directional one: z is in the direction of the sun's motion, x is in the radial 
direction and x is perpendicular to the axis of the galaxy. The asymmetry is also given. 
Qmin = was assumed. 



isothermal caustic rings 



A dir. 
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asym + 


asym 


z 


0.02 0.50 


0.92 0.75 0.25 


0.50 
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0.16 0.16 


0.22 


1.00 
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0.16 0.16 


0.37 0.24 


0.21 


1 z 


0.04 0.58 


0.90 




1 y 


0.12 0.12 







1 X 


0.17 0.17 








Fig. |^ and Fig. ||. We must emphasize that the two graphs of Fig. || do not 
contain the entire dependence on the LSP mass. This is due to the fact that 
there is the extra factor rnZ 1 in Eq. (FL9) and a factor of fj% arising from Si, 
i = S, spin, V, see Eqs (Q), (0), and (|10[). All these factors combined lead to 
a constant. There remains, however, an LSP mass dependence, which is due to 
the fact that the nucleon cross section itself depends on the LSP mass. 

Figs I and I were obtained for the scalar interaction, but similar behavior 
is expected for the spin contribution. The scale of the total spin contribution, 
however, is going to be very different. We should also mention that in the non 
directional experiments the modulation 2h\ is small, .i.e. for A = less than 4% 
for Qmin = and 12% for Qmin = 10 KeV (at the expense of the total number 
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of counts). For A = 1 there in no change for Q m in — 0, but it can go as high as 
24% for Q min = 10 KeV. 

For the directional rates It is instructive to examine the reduction factors 
along the three axes, i.e along +z, —z, +y, —y, +x and —x. These depend on the 
nuclear parameters, the reduced mass, the energy cutoff Q m in and A ]27| |. Since 
f re d is the ratio of two parameters, its dependence on Q m in and the LSP mass 
is mild. So we present results for Q m i n — and give their average as a function 
of the LSP mass (see Table |l|). As expected the maximum rate is along the sun's 
direction of motion, i.e opposite to its velocity (—z) in the Gaussian distribution 
and +z in the case of caustic rings. In fact we find that n(— z) is around 0.5 (no 
asymmetry) and around 0.6 (maximum asymmetry, A = 1.0). It is not very dif- 
ferent from the naively expected f rec i = 1/ (27r) = n = 1 . The asymmetry along 
the sun's direction of motion, asym = \Rdir{—) — Rdir{+)\/ {Rdir{~) + Rdir{+)) 
is quite characteristic and, for a given direction, it depends on the velocity dis- 
tribution. It is large in the direction of the sun's motion in the isothermal models 
and a bit smaller in caustic rings. The rate in the other directions is quite a bit 
smaller (see Table 0). In this case for the isothermal model the asymmetry is 
zero, while in the case of caustic rings it is sizable. So the two models can be eas- 
ily distinguished. The disadvantage of smaller rates in the plane perpendicular 
to the sun's velocity may be compensated by the very large and characteristic 
modulation. It is interesting to note that even the functional dependence on the 
phase of the Earth a changes substantially with the direction of observation. 

In conclusion: in the case of directional non modulated rates we expect un- 
ambiguous correlation with the motion on the sun, which can be explored by 
the experimentalists. The reduction factor in the direction of the motion of the 
sun is approximately only l/(47r) relative to the non directional experiments. In 
the plane perpendicular to the motion of the sun we expect interesting modula- 
tion signals, but the reduction factor becomes worse. These difficulties may be 
reduced, if the TPC counters can make observations in many directions simul- 
taneously. 

This work was supported by the European Union under the contracts RTN 
No HPRN-CT-2000-00148 and TMR No. ERBFMRX-CT96-0090. 
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